  node_number=100;
 netarea=10;communicationrange=2.5;
 interface=3;
%  channel_number=30;

channelpool=10:2:30;
 I=zeros(1,length(channelpool));D=zeros(1,length(channelpool));H=zeros(1,length(channelpool));RUNTIME=zeros(1,length(channelpool));
rank=0;
 for channel_number=channelpool
    rank=rank+1;
for i=1:10
    [runtimes,In,Dn,Hn]=DM_DTMalg(netarea,communicationrange,node_number,interface,channel_number);
    I(rank)=I(rank)+In;D(rank)=D(rank)+Dn;H(rank)=H(rank)+Hn;RUNTIME(rank)=RUNTIME(rank)+runtimes;
end
end
I/10
D/10
H/10
RUNTIME/10
plot(channelpool,I/20,'r')
plot(channelpool,D/20,'g')
plot(channelpool,H/20,'b')
plot(channelpool,RUNTIME/20,'y')
